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Abstract 

The variation with microwave frequency and temperature of previously reported anomalous peaks 
in the EPR spectra of Mni2-acetate, under large transverse helds, reveals that the molecular easy 
magnetization axes are tilted with respect to the global symmetry direction. More importantly, on 
the basis of the angle-dependence of fine structures observed in the EPR spectra we infer that the 
tilt distribution must be discrete, as was previously suspected from studies which demonstrated 
the presence of a locally varying rhombic anisotropy [S. Hill et al, Phys. Rev. Lett. 90, 217204 
(2003)]. The tilts are confined to two orthogonal planes, and the distribution extends up to ~ 1.7° 
degrees away from the the global easy (z-) axis. We ascribe the tilting to the hydrogen-bonding 
effect associated with the disordered acetic acid solvent molecules. The effect is considerably larger 
than deduced from x-ray diffraction analyses. These data constitute the sought-after evidence for 
the presence of transverse fields in Mni2-acetate, and provide a possible explanation for the lack 
of selection rules in the resonant quantum tunneling behavior seen in low-temperature hysteresis 
experiments for this S = 10 system. 

PACS numbers: 75.50.Xx, 75.60.Jk, 75.75.+a, 76.30.-v 
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I. INTRODUCTION 



Since the discovery of magnetic quantum tunneling (MQT) in 
[Mn 12 Oi 2 (CH 3 COO)i 6 (H 2 0) 4 ]- 2CH 3 COOH-4H 2 (Mn 12 -ac),i2^^ single molecule 
magnets (SMMs) have become the focus of considerable experimental and theoretical 
interest due to their novel quantum properties^ and possible future use in quantum 
computational devices.—^ Mn 12 -ac remains the most widely studied SMM due to its 
large spin ground state (S = 10, see refs ), together with a considerable easy-axis 

magneto-crystalline anisotropy.i& These combined factors result in a sizeable kinetic barrier 
against spin reversal at the molecular level, leading to slow magnetization relaxation and 
hysteresis (bistability) at low temperatures (below ~3 K)££ When a DC magnetic field 
is applied parallel to the easy-axis of a single crystal of Mni 2 -ac, sharp steps are observed 
in its hysteresis loops at well defined field strengths.— i 6 i 7 i 8 i n The enhanced magnetization 
relaxation at these steps is the result of resonant MQT. 

While a clearer picture concerning the mechanism of MQT in Mni 2 -ac is beginning to 
emerge , 12 i 13 i 14 i 15 many interesting problems remain. In particular, current theoretical models 
assume the presence of quadratic and quartic transverse crystal-field interactions in the 
spin Hamiltonian [0| = + S?_) and Of = 1(5"+ + S*)], where the former has been 

ascribed to solvent disorder.- However, these interactions, which contain only even powers 
of the raising and lowering operators, do not provide an explanation for the observation 
of 'odd' MQT steps in the hysteresis loops, i.e. tunneling via resonances between levels 
whose spin projections (m s ) differ by an odd integer.^ It is generally recognized that the 
underlying mechanism must involve internal transverse fields , 12 ' 1 ^ but their source has yet 
to be identified. Prokof 'ev and Stamp 17 proposed that hyperfine interactions might provide 
the answer; however, the calculated MQT rates were too small. Subsequently, Garanin and 
Chudnovsky 12 suggested that strains or structural defects in the crystal lattice (dislocations) 
could lead to a distribution of tilts of the magnetic easy-axes at the Mn i2 -ac cluster sites; 
upon application of an external magnetic field, such tilts would result in a distribution of 
transverse fields, even when the field is applied parallel to the global symmetry (z-) axis of 
the crystal. While this model was in qualitative agreement with the fact that experiments 
revealed significant distributions in the finite field MQT rates , 18 > 19 the measured distributions 
were considerably narrower than those predicted by Garanin and Chudnovsky.— 
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Although some preliminary spectroscopic data have been reported in support of 
dislo cat ions, ~ a more plausible model leading to easy-axis-tilting was recently suggested 
by Cornia et al.— Based on a detailed analysis of existing x-ray diffraction data and some 
approximate electronic structure calculations, they propose that such a local symmetry 
lowering and easy-axis-tilting can be ascribed to the presence of the two acetic acids of crys- 
tallization in the unit cell. Experiments in support of this model were recently published 
by us, 14 and by del Barco et al.— However, evidence for the easy- axis-tilting was lacking 
in these investigations, though we note that tilts have recently been reported for a related 
Mni2 complex in a separate study by del Barco et al.— Here we report electron paramagnetic 
resonance (EPR) measurements on precisely (in-situ) aligned Mn 12 -ac single crystals under 
large applied fields in the transverse direction. At relatively low frequencies (< 90 GHz), 
as the field is rotated away from the hard plane, simulations show that the EPR intensity 
should oscillate between two series of resonances excited from "even" and "odd" spin states 
(labeled a and (3 respectively, see following section). The magnetic dipole matrix elements 
and transition frequencies associated with these "even" and "odd" resonances are extremely 
sensitive to the field orientation (providing < 1° resolution), allowing for a precise appraisal 
of the Cornia model. Indeed, experiments reveal a significant overlap of the a and /3 res- 
onances, providing conclusive evidence for the presence of tilts. However, we find that the 
magnitude of the tilts is a factor of 4 — 6 larger than predicted by Cornia et al.,— which 
should be of both theoretical and experimental significance in regards to understanding the 
spin dynamics of Mni2-ac. Our conclusions also provide an explanation for several previ- 
ously reported anomalous EPR transitions (labeled /^-resonances) which cannot be explained 
within the widely accepted giant spin (S = 10) model described below. 22 

We emphasize from the outset that our conclusions are based upon the angle-dependent 
behavior of EPR fine structures which have previously been explained^ 4 - within the context 
of the Cornia solvent-disorder picture,^ and have been independently supported by magne- 
tization measurements.- 1 ^ The easy axis tilting is inferred from the persistence of EPR peaks 
for field orientations which are tilted significantly away from the angles where the intensity 
should have vanished if all molecules were aligned. The discrete nature of the tilting can then 
be seen from the fact that the different fine structures exhibit distinct angle dependencies. 
This behavior is completely reproducible in Mni2-ac samples prepared by different methods 
and by different groups . 2 ?' 24 Furthermore, our method of analysis is not sensitive to the 
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intricacies associated with lineshape analyses, which could be influenced by complex many- 
body effects, e.g. magnetic spin-spin interactions.— Finally, we note that all of the unusual 
solvent-induced anomalies reported in this paper for Mn 12 -ac are absent in the EPR spectra 
for a related high-symmetry (5 4 ) Mn 12 complex [Mn 12 12 (0 2 CCH 2 Br) 16 (H 2 0) 4 4CH 2 Cl 2 ] 
which possesses a full compliment of four CH 2 C1 2 solvent molecules per Mn 12 .— 

II. BACKGROUND 

The basis of our experiment is illustrated in Fig. 1, which shows the spin energy level 
diagram for Mni 2 -ac (Fig. la), based on an exact diagonalization of the standard S = 10 
spin Hamiltonia n 8 i 13 i 14 i 15 i 21 i 22 i 27 i 28 i 29 under a transverse magnetic field B, and assuming £4 
symmetry: 

H = DS 2 z +iJ, B B-g-S+B26 < l + Bidj, (1) 

Here, S z is the projection of the spin operator S along z, and D (< 0) is the uniaxial 
anisotropy parameter; the second term is the Zeeman interaction; and the remaining terms 
represent higher order crystal field interactions (the operators 0° and 0\ have their usual 
meaning).^ 8 . We have used acceptable parameters for the simulations in Fig. 1: D — —0.455 
cm -1 , B® = —2.0 x 10 -5 cm -1 and B\ = ±3.2 x 10~ 5 cm -1 . We note that this parameter 
set yields the optimum agreement with a large body of single crystal EPR data , 14 ' 30 ' 31 ' 32 
including the present study, which additionally considers the effects of .E-strain and easy axis 
tilting. At this stage, we do not explicitly include lower symmetry interactions in Eq. 1 due to 
disorder, e.g. second-order rhombicity, or easy-axis tilting. We begin by considering a single 
molecule and the effect of field-mis-alignment away from the hard plane. We then compare 
experimental data with simulations in order to quantify the easy-axis tilting caused by the 
solvent-disorder in Mn 12 -ac. Finally, we consider the rhombic anisotropy associated with this 
disorder. We note that an account of the influence of the solvent-disorder-induced rhombic 
anisotropy [E(S^. — 5* 2 )], and its affect on the transverse field EPR fine structure, has been 
presented previously by us . 14 ' 31 In the following, we use polar coordinates to parameterize 
the field orientation: 6 represents the angle between the applied field and the global easy 
(z-) axis of the single crystal; cj) represents the azimuthal angle, i.e. the angle between the 
intrinsic hard four-fold (x-) B\ axis and the projection of the applied field onto the hard 
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plane. 

In the high-field limit (g/^eB > l-C^S), and for the standard EPR geometry (microwave 
field Hi -LB), one expects a total of 20 Am = ±1 EPR transitions between the 21 (= 2S + 1) 
spin-states for S = 10. For clarity, only a few of the lowest-lying levels are shown in Fig. la 
for the case of a magnetic field applied precisely along the medium four-fold axis of a single 
molecule (i.e. 8 = 90° and = 45°); these levels are labeled on the right-hand-side of Fig. la 
according to an S x basis, where m x (= integer) is the projection of the total spin along the 
applied field axis. The high-field EPR spectra are dominated by transitions between adjacent 
levels, i.e. Am x = ±1. If one follows the EPR spectra to lower fields, into a region where 
the Zeeman and axial terms in Eq. 1 become comparable (g/i^B ~ |-D|S0, one finds that the 
transitions may be grouped into two categories: i) those between levels which evolve from 
m z = ±i (i = integer) zero-field doublets, which we label a (represented by blue sticks in 
Fig. la); and ii) those between levels which evolve from adjacent zero-field doublets, which 
we label j3 (represented by red solid circles in Fig. la). This distinction between a and j3 
resonances is based on a zero-field S z basis, where m z is the projection of the total spin along 
the uniaxial direction of the crystal — the levels are labeled according to this convention 
on the left-hand-side of Fig. la. Throughout the remainder of the paper, we number all 
transitions according to the absolute value of m x (high-field S x representation) associated 
with the level from which the transition was excited, preceded by either a or (3 (low-field S z 
representation) in order to distinguish between the two categories of resonances. Therefore, 
the highest field blue stick in Fig. la corresponds to alO, while the highest field red circle 
corresponds to j39. 33 

In order to set the scene, we first review the status of earlier single-crystal EPR studies. 
The observation most pertinent to the present work is seen in Fig. lb, which compares various 
calculations of the a and (3 transition frequencies with actual 8 K EPR peak positions . 30 ! 31 
For these earlier experiments, data were obtained with the field applied parallel to the hard 
plane of a single-crystal sample to within an accuracy of ~ 1.5° (6 = 90° ± 1.5°). Within 
the hard plane, a crude attempt was made to align the magnetic field along the medium 
four-fold axis. Based on our more recent angle dependent studies on accurately aligned 
single crystals (section III and ref. Q|), we estimate that such an alignment was achieved 
to within an accuracy of approximately 10° (i.e. <p ~ 45° ± 10°). For the purposes of 
the following discussion, any mis-alignment within the hard plane is not important. In 
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the absence of tilting, the calculated a transition frequencies decrease smoothly to zero 
(thick blue curves in Fig. lb), whereas the (3 transition frequencies (thick red curves in 
Fig. lb) go through minima which are on the order of 90 GHz for the three highest field 
branches. The actual /^-transition data, on the other hand, deviate significantly from these 
predictions, i.e. they do not exhibit a minimum frequency, but instead follow monotonic 
curves to the lowest frequencies investigated. The same trend has been noted for field 
alignment along different directions within the hard plane . 30 ' 31 ' 32 For comparison, Fig. 2 
displays representative experimental EPR spectra obtained with the field approximately 
along the hard four-fold axis (~ cb = 0°), also within ±1.5° of the hard/medium plane 
(6 = 90° ± 1.5°);— note that the /3-resonances are again seen to the lowest frequencies 
studied (44 GHz). 

Following our more recent angle-dependent studies that provided clear indications for a 
breakdown of the four-fold symmetry of Mn 12 -ac,— and in light of the Cornia model,— we 
conjecture that the anomalous behavior of the f3 transitions might be explained in terms 
of tilts of the easy-axes of magnetization at a local scale. To better illustrate this idea, 
we have included tilted field calculations in Fig. lb (thin curves: blue— a, red—/?) for a 
few of the highest field transitions. Each tilted field curve belongs to a family of curves 
which deviate successively (in 0.5° increments) from one of the zero-tilt curves (9 = 90°, 
thick lines); the tilted field curves have been truncated at low fields in order to simplify the 
figure. Tilting the field away from the hard plane results in a field component parallel to 
the easy axis. This "longitudinal" field (B\\) has the effect of lifting the degeneracy of the 
of the m z = ±i zero-field doublets in zeroth order. 34 Consequently, instead of tending to 
zero-frequency, the a resonance frequencies tend to successively larger offsets as the field 
is tilted away from the hard plane, as seen by the thin blue curves and indicated by the 
blue arrow (for alO) in Fig. lb. The opposite is true for the f3 transitions. Since the 
longitudinal field splits the low-field m z = ±z doublets, this results in a reduction in the 
closest approach of levels in adjacent doublets — hence, to a reduction in the minimum 
frequency of the (3 transitions; again, this behavior is born out by the thin red curves and 
indicated by the red arrow (for /39) in Fig. lb. Thus, tilting could explain the continuation 
of the (3 resonance data to frequencies well below the theoretical minima given by the thick 
red curves (8 = 90°, or B» = 0) in Fig. lb. However, the observed trend cannot be ascribed 
to a simple mis-alignment of the sample, because the alO (highest field blue squares) and 
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(39 (highest field red circles) resonances are both seen to frequencies below 50 GHz (see also 
Fig. 2). The observation of /39 at 45 GHz would imply a sample mis-alignment of 2 — 3°, 
which is completely incompatible with the observed behavior of the alO data which tracks 
the thick blue curve (9 = 90°, or B» = 0) in Fig. lb; note that each successive thin curve 
corresponds to an additional 0.5° of tilt away from the hard plane. Thus, it would appear 
that the data in Figs. 1 and 2 reflect several molecular orientations. 

A rigorous comparison between experiment and Eq. 1, in this frequency/field range where 
the Zeeman and axial terms are comparable, requires calculation of magnetic dipole matrix 
elements, i.e. a full simulation of the EPR spectra. Figs. 3 and 4 show 15 K EPR simulations, 
as a function of the polar angle 9 between the applied field and the easy axis of a single 
molecule, and for two of the lowest frequencies used in our experiments (44 GHz and 62 GHz); 
the simulations are limited to the 6 = 80° to 90° range in 0.2° increments. In Figs 3a and 4a, 
the vertical scale represents microwave absorption while, in Figs 3b and 4b, absorption is 
indicated by the darker shaded regions. The simulations were generated using the software 
package SIM; ?^ no rhombic (E) term was included at this stage, and the simulations 
assume perfect alignment of all molecules in the sample, i.e. no orientational averaging was 
employed. As expected, the simulations confirm that only the a peaks are observed for fields 
parallel to the hard plane (6 = 90°), for / < 90 GHz. However, tilting the field away from the 
hard plane leads to a very abrupt suppression of the magnetic dipole matrix elements (EPR 
selection rules) for the a transitions — this behavior was not obvious from Fig. 1. In just 1° 
of rotation, alO vanishes completely at 44 GHz (~ 1.5° for 62 GHz). This is then followed 
by a range of almost 2° where neither alO nor (39 are observed at the lowest frequency of 
44 GHz, i.e. (39 does not appear until the field is tilted about 2.8° away from the hard plane. 
Similar trends are seen at both frequencies for a8 and (37, etc., albeit over differing angle 
ranges. Remarkably, the simulations are suggestive of a symmetry effect in the transverse- 
field (-Bj_) EPR spectra, as a function of the longitudinal field component (Si < B±). 
The dashed curves on the contour plots in Figs. 3b and 4b indicate lines of constant B\\, 
demonstrating that all of the a resonances disappear (and all of the (3 resonances appear) at 
roughly the same longitudinal field strengths. Subsequently, the (3 resonances disappear and 
the a' resonances re-appear. This behavior is reminiscent of the tunnel-splittings seen in the 
Fes SMM as a function of transverse field, for different static longitudinal fields.— Indeed, 
closer inspection of Fig. lb reveals that the (3 transition frequencies tend to zero at certain 
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values of B\\, while the simulations in Figs. 3 and 4 suggest that the (3 transitions matrix 
elements tend to zero as well at these same angles. Thus, the analogy with the quenching 
of the ground-state tunnel splittings in Fe§ is significant, albeit the degeneracies involve 
excited levels evolving from different zero-field doublets in the present case. We shall not 
pursue this analogy further here. However, we note that the (3 series of resonances is exactly 
what one would expect from a spin S = 9 system in a perpendicular field with the same 
Hamiltonian parameters as the parent S — 10 system ; 26 i 30 i 38 meanwhile a8' corresponds to 
the first expected peak in a series belonging to a spin S = 8 system, and so on. Furthermore, 
all resonances except for alO occur from excited levels. Thus, the tilted field data have the 
appearance of spectra from excited multiplets having successively lower total spin, as was 
originally proposed by us as an explanation for the anomalous (3 resonances. 30 

An intensity analysis of the measured low-frequency (< 90 GHz) (39 peaks in Fig. 2 yields 
an (approximate) activation energy consistent with an excitation from the m x = —9 state 
within the S — 10 manifold. We should point out that, in order to precisely fit the (39 
intensity, one must first fully comprehend its origin (we save such an analysis until the end 
of this article). The clearest indication that f39 likely occurs within the S = 10 manifold 
can be seen by noting that the temperature dependencies of the (39 and a8 transitions in 
Fig. 2 are very similar, suggesting that they both involve energy levels which are close in 
energy. Thus, in the following sections, we do not consider the possibility of excitations 
within excited state (S ^ 10) manifolds. However, the simultaneous observation of alO and 
(39 down to a frequency of 44 GHz (Fig. 2) is not consistent with the simulations in Fig. 3, 
which assume aligned molecules. For this reason, we argue that the molecular easy-axes are 
not perfectly aligned, i.e. we propose a distribution of easy-axis alignments centered about 
the global four-fold axis of the crystal. This would result in overlapping angle dependent 
features (such as those in Figs. 3 and 4) from different parts of the distribution, i.e. for 
a given field orientation one may observe (3 peaks from tilted (aligned) molecules together 
with a peaks due to aligned (tilted) molecules. As already discussed, the origin for such 
a distribution can be understood in terms of solvent disorder. Indeed, the calculations by 
Cornia et al. predict easy-axis tilts of up to 0.5°.— The absence of either alO or (39 over an 
approximately 2° range in the simulation in Fig. 3 suggests that the tilt distribution extends 
at least ±1°, since both alO and (39 are observed in the 44 GHz spectra (Fig. 2). 
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III. EXPERIMENTAL 



In order to confirm the above hypothesis, we have recently carried out extremely pre- 
cise angle-dependent measurements for 9 rotations away from the hard plane (~ 0.1°). In 
contrast to earlier angle- dependent investigations within the hard plane (</> rotations),— 
which were achieved using a rotating split-pair magnet, the present study required fields 
exceeding those available in the split-pair. Consequently, a unique cavity was designed and 
constructed, enabling in-situ rotation of the sample on the end-plate of the cylindrical cavity 
with an angle step of 0.18°, and in fields of up to 45 tesla (at the National High Magnetic 
Field Laboratory). Details concerning this cavity will be published elsewhere,— and an ac- 
count of our EPR spectrometer is published in ref. The present investigation was carried 
out in a standard 9 tesla superconducting solenoid at the University of Florida. The sample 
was positioned on the cavity end-plate so that the plane of rotation coincided approximately 
with one of the four large flat faces of the needle-shaped sample. Based on comparisons with 
magnetic measurements— we believe that this plane of rotation is inclined approximately 
10° with respect to a plane containing the easy-axis and one of the hard four-fold axes, i.e. 
we believe cb ~ —10° in these investigations (see Fig. 6 below). 32 Single crystal samples of 
[Mn 12 Oi 2 (CH 3 COO)i 6 (H 2 0)4]- 2CH 3 COOH-4H 2 were grown using literature methods.i 3 - 
Magnetic field dependent absorption was recorded at 15 K and at a frequency of 61.9 GHz, 
for field orientations roughly 3.5° either side of the hard plane (9 = 86.5° — > 93.5°), in 0.18° 
increments. The experimental spectra displayed in Fig. 5 are presented in several different 
ways in order to aid direct comparison with Fig. 4. We note that the restricted bore size of 
our high-field magnet forced a lower bound on the resonance frequencies possible in the spe- 
cially designed rotating cylindrical cavity— Thus, high-field measurements below 61.9 GHz 
were not possible. 

Immediately apparent from Fig. 5 is the fact that the a and (3 peaks overlap over a 
substantial angle range (~ 2° — 3°) — particularly culO, (39, and a.8 — whereas this is not the 
case in Fig. 4, thus providing direct confirmation for the tilt distribution. One also observes 
a weak (39 peak in Fig. 5b for field alignment precisely within the hard plane (9 = 90°); in 
fact, (39 is reduced below the 15% level for only about 0.5° either side of the hard plane, in 
comparison to its complete absence over a 90° ±2° range in Fig. 4b. All of these facts suggest 
that the easy axes of some of the molecules must be tilted by up to the order of ±1.5°, i.e. 
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substantially greater than the tilts predicted by Cornia et al.— Either that, or there exists a 
significant distribution of internal transverse fields within the sample. However, as previously 
discussed, this has been ruled out by previous authors.— Indeed, internal fields are rather 
weak , 40 i 41 and cannot explain the present observations. Even for the worst possible sample 
shape, depolarization effects could account for no more than a 0.2° distribution of field 
orientations at the applied field strengths employed in these investigations. A more subtle, 
yet significant, aspect of the data in Fig. 5 is the presence of shoulders on the high- field sides 
of most of the EPR peaks (see Fig. 5a). These shoulders have been emphasized in Fig. 5c 
via the 50% contour lines (this corresponds to 50% of the maximum absorption). Such 
shoulders have been noted previously by us,— and ascribed to the solvent-disorder-induced 
rhombic anisotropy first proposed by Cornia et al.— It is evident from Figs. 5a and c that 
the shoulder persists over a narrower angle range, when compared to the main part of the 
EPR absorption peaks, i.e. suppression of the shoulder relative to the main peak is quite 
abrupt, occurring close to 9 = 90 ± 1° for alO. The occurrence of distinct regions in the 3D 
absorption plot (Fig. 5c), with distinct angle dependencies, is suggestive of a discrete form 
of disorder, as originally suggested in our earlier EPR investigations.— As will become 
apparent later in this article, the shoulder is attributable to a finite rhombic term. However, 
the angle dependence cannot be ascribed to a rhombic term alone — one must additionally 
consider the combined effects of easy-axis tilting and rhombicity. In doing this, one may 
understand the distinct angle dependence of the main peaks and the shoulders (Fig. 5) as 
being due to molecules tilted in distinct (orthogonal) planes. 



IV. DISCUSSION 



The origin of two separate angle dependencies finds a natural explanation if one assumes 
that the easy-axis tilting and rhombicity are connected (as predicted by Cornia), and that 
the tilting is confined to directions determined by the principal axes of the rhombic zero- 
field tensor. We now assume that two of the Mn 12 variants originally discussed by Cornia 
et al. (n = 1 and n = 3)^ have their easy axes tilted significantly, and that these tilts are 
confined to two orthogonal planes defined by the associated rhombic zero-field tensor and 
the global four-fold (z-) axis. The n = 1 and n = 3 variants comprise 50% of the total 
number of molecules in Cornia's model, and both species are predicted to have significant 
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E values (1.6 x 1(T 3 cm- 1 ), 13 Of this 50%, half of the molecules (25% of the total - the 
A molecules) are expected to have their hard two-fold axes aligned along a single direction 
within the hard plane. Then, because there exist four equivalent positions for the acetic acid 
of crystallization related by a four-fold rotation about z, the remaining half of the n = 1 
and n = 3 species (25% of the total — the B molecules) will have their hard axes aligned 
90° away from the hard axes of the A molecules, i.e. parallel to the medium two-fold axes 
of the A molecules. This situation is depicted by the polar (</>-) plot in the lower panel of 
Fig. 6. The two-fold hard axes are denoted HE, with the subscript A or B to distinguish 
between the two sub-sets of molecule. The orientations of these axes relative to the four-fold 
hard-axes [HBf at = 0°, 90°, 180°, and 270°) is based on previous angle- dependent EPR 
investigations for rotations within the hard plane, where it was determined that the principal 
axes of the E and B\ tensors were mis- aligned by 30° . 14 i 32 Next, we assume that the easy axis 
tilts are also confined along these two distinct /discrete directions, as indicated by the shaded 
regions in the lower panel of Fig. 6. In other words, the well defined rhombic distortion and 
molecular easy axis tilting are directly related. Based on the ensuing analysis, we find that 
molecules tilted in one of the two orthogonal planes have their HE axes perpendicular to 
that plane, as shown in the upper panel of Fig. 6. Thus, the A (B) molecules are tilted in 
a plane containing the HEb (HEa) axis. 

It is now possible to see how one might expect two distinct angle dependencies for different 
parts of the EPR spectrum, i.e. the shoulders and the main peaks. The thick purple arrow 
along (ft = —10° (and 170°) in the lower panel of Fig. 6 represents the best estimate (vide 
infra) of the plane of rotation for the EPR spectra presented in Fig. 5. This plane is inclined 
closer to HE A (0 = -30° and 150°) than to HE B (0 = 60° and 240°). Therefore, for fields 
close to the hard plane (9 = 90°), the rhombic term (EO2) associated with the n = 1 and 
n = 3 variants will lead to different EPR peak-position-shifts for the A and B molecules, 
thereby enabling one to resolve their contributions to the EPR spectra for a sufficiently 
large E value. In our experiments, the plane of rotation is closer to the tilt-plane for the B 
molecules than for the A molecules (A</> = 20° as opposed to 70°). Thus, the projection of 
the tilt distribution onto the plane of rotation is correspondingly greater for the B molecules 
than for the A molecules. The magnitudes of these projections scale as the cosine of the 
angle (A</>) between the field rotation plane and the corresponding tilt-plane. Indeed, as 
explained later, it is on this basis that we deduced the orientation of the plane of rotation 
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relative to HE a and HEb for the experiments presented in Fig. 5. The molecules having 
the greater projection of tilts onto the rotation plane will, thus, be expected to exhibit the 
greater spread in their associated EPR spectra away from the hard plane. In this scenario, 
one may attribute the shoulders on the high-field sides of the main peaks in Fig. 5 to the 
A molecules in Fig. 6, since their tilt distribution will have a narrower projection onto the 
rotation plane (see projections onto field rotation plane represented by blue and red dashed 
lines in Fig. 6). Furthermore, the fact that the rotation plane is close to the hard axes of the 
A molecules is fully consistent with the appearance of the shoulder on the high-field sides 
of the main peaks (assuming E > 0). Meanwhile, the contribution of the B molecules is 
buried within the main peak, as explained below. 

Based on the geometry depicted in Fig. 6, we set out to simulate the data in Fig. 5. The 
discussion above concerns only 50% of the possible variants in Cornia's model,— i.e. the two 
dominant species possessing a significant i?-term. Furthermore, only half of these (25% of 
the total) are tilted in the plane of rotation. For the sake of simplicity, we assume that the 
remaining 50% of the molecules (C) do not possess a significant lvalue. This is a slight 
simplification of Cornia's model, since it is known that the n = 2 trans species (comprising 
12.5% of the molecules) are predicted to have an .E-value comparable to the n — 1 and n = 3 
variants. 13 However, one should first recall that Cornia's calculations are approximate, and 
that our previous EPR studies have indicated that the magnitudes of the .E-values may 
have been underestimated.— Second, introducing too many different species would result in 
an over-parameterization of our model. The main purpose of the following simulation is to 
demonstrate that one can explain essentially all aspects of the single-crystal EPR spectra 
for Mn 12 -ac based on a simple model involving three just three molecular species, together 
with discrete easy-axis tilting, in the spirit of the original Cornia proposal. It is not our 
intent that this model be viewed as definitive. Fig. 7 shows three gray-scale contour plots 
representing the expected angle (9-) dependent contributions to the EPR intensity for the 
three species A, B, and C; they assume the same experimental conditions as Figs. 4 and 
5, i.e. T = 15 K and a frequency of 62 GHz. The simulations in Figs. 7a and b assume 
that the A and B molecules are tilted respectively along the HEa and HEb directions 
(0 = —30° and 60°) i.e. the orientations of the tilt planes are "discrete." Meanwhile, we 
assume a non-discrete distribution of the tilt angles along these two directions, which cuts-off 
at 9 = 90° ±1.7°, as indicated by the shaded region in the lower panel of Fig. 6. In addition, 
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we assume an .E-value of 0.008 cm -1 for the A and B species. The simulation in Fig. 7c 
assumes no _E-value. Consequently, the EPR intensity due to the C molecules contributes 
to the bulk of the central portion of the summed EPR peaks (Fig. 8). We included a small 
random (i.e. not confined to planes) distribution of tilts (±1°) for the C molecules as a 
means of taking into account possible variations in the zero-field parameters/orientations 
associated with the four Cornia variants (n = 0, n = 4, n = 2 cis and trans) which comprise 
this remaining 50% of the molecules. It is the finite lvalue associated with the A and 
B species which enables us to resolve the discrete nature of their tilting directions. We 
cannot rule out similar properties for the remaining species (C), but their smaller E- values 
do not allow us to observe such discrete behavior. Thus, the choice of random tilts for 
the C molecules is quite arbitrary. Nevertheless, this distribution provides better overall 
qualitative agreement with the observed out-of-plane angle dependence in Fig. 5. 

Agreement between experiment (Fig. 5) and the subsequent summation of the simulated 
spectra (Fig. 8) is rather sensitive to the precise cut-off angle of 6 = 90° ± 1.7° for the tilt 
distribution of the A and B molecules; however, it is not so sensitive to the exact shape 
of the distribution. For example, a discrete distribution with the A and B molecules all 
having precisely 9 = 90° ±1.7° tilts along HE a and HEb reproduces many aspects of 
the experiments. However, such simulations contain many more fine structures which are 
averaged out by considering a smoother distribution such as the one employed in the Fig. 8 
simulations. Once again, we stress that this parameter set should not be viewed as definitive. 
However, we emphasize that, in order to simulate the main qualitative trends in the data 
(e.g. the distinct angle dependencies of the main peaks and high- field shoulders), a discrete 
tilt distribution is necessary. As discussed above, the simulations are not so sensitive to 
the shape of the tilt distribution within each tilt plane. Thus, it is not easy to make direct 
comparisons with the distributions of tilts inferred by other methods.— However, we note 
that only a small subset of the molecules (25%), corresponding to half of the lowest symmetry 
Cornia variants^ are able to contribute to the anomalous EPR intensity in our experiment. 
Of these molecules, we can conclude that fewer than half (~ 10% of the total) are tilted 
by more than 1°. While comparisons between the gray-scale plots are not so sensitive to 
the .E-value, comparisons between simulated and actual experimental spectra are extremely 
sensitive to E (see Fig. 10 below); the value of 0.008(2) cm -1 gives the best agreement with 
a wide body of single- crystal data collected over several years . 14 i 22 i? ( - ) i? 1 i 44 This E value is also 
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in excellent agreement with a recent first principles calculation by Park et al.— 

As discussed above, the relative angle (A0) between the field rotation plane and the A and 
B tilt-planes (HEb and HE a respectively) was chosen so that the EPR peaks corresponding 
to the A and B molecules (Figs. 7a and b) cover the same angle ranges as the high-field 
shoulders [9 ~ 90° ± 1° for alO) and the central portions of the peaks (9 ~ 90° ± 2° for 
alO) respectively in Fig. 5. The best agreement was obtained with the field rotation plane 
oriented A0 = 20° away from HE a (Fig. 6). Thus, we estimate that the HE a and HEb 
axes are oriented =f20° and ±70° with respect to the sample faces, while the hard four-fold 
axes are roughly ±10° away (see Fig. 6), which also agrees with magnetic and structural 
data . 15 i 23 The simulated spectra obtained according to the aforementioned procedure were 
then normalized and summed in the ratio A:B:C = 25% : 25% : 50%. The results of this 
summation are displayed in Fig. 8. Agreement with the experimental data in Fig. 5 is quite 
impressive. First and foremost, the simulations including tilts account for essentially all of 
the anomalous aspects of the data which could not be explained without tilts (Fig. 4). For 
example: the alO and (39 peaks overlap by about 2° between roughly 9 = 90° ± 0.5° and 
90° ± 2.5°; the (39 peak disappears (below the 10% level) over a roughly ±0.5° range either 
side of 9 = 90°. There is also clear evidence for a shoulder on the high-field sides of the main 
EPR peaks, although this shoulder is more apparent in Fig. 5a than 5b. The overall widths 
and shapes of the resonances are also reproduced fairly well, including the broad low-field 
tail and the fairly abrupt decrease in intensity on the high-field sides of the peaks (a 0.1 tesla 
Gaussian lineshape was employed in the individual simulations in Fig. 7). We note from 
Fig. 5 that an extremely weak (39 signal persists even for fields parallel to the hard plane. 
We were unable to reproduce this behavior in the simulations while, at the same time, 
maintaining good agreement with other aspects of the experimental data, e.g. the angle 
range of overlap between the a and (3 resonances. We therefore speculate that this weak 
remnant signal may be due to excitations within a higher lying S = 9 state, as previously 
proposed by usM However, more careful temperature dependent studies will be required 
in order to establish the location of this excited state relative to the S = 10 ground state. 
Indeed, recent measurements on a related high symmetry Mn 12 species (without significant 
easy axis tilting) have enabled precisely such an analysis. 26 

Having satisfactorily reproduced the spectra for out-of-plane rotations (Figs. 5 and 8), 
we consider previously published spectra obtained as a function of the field orientation 
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within the hard plane.— One of the puzzling questions concerning these earlier measurements 
concerned the observation of only a single satellite peak (shoulder) on the high-field side of 
resonances for some field orientations. The modulation of the widths of the EPR peaks was 
originally explained in terms of the Cornia model.— So why does one never see both low-field 
and high-field shoulders, corresponding to the positive and negative EPR line-position-shifts 
induced by the different signs of the rhombic term? We note that, for a field along HEa 
(HEb), one expects an upward (downward) shift in field of the absorption due to the A 
molecules, together with a downward (upward) shift due to the B molecules. The apparent 
explanation is illustrated in Fig. 9, which displays the evolution of absorption due to one of 
the tilted species (either A or B) as a function of the field orientation (within the hard plane) 
relative to the hard two-fold axis for that species; note — this is not the same as the angle cf>. 
With the field along the hard two-fold axis HE (<pE = 0° in Fig. 9), the resonance occurs at 
the highest field position, while the projected distribution of tilts onto the rotation plane is 
at its narrowest. Consequently, the EPR peak is narrow. Conversely, with the field along the 
medium two- fold axis ((f) e = 90° in Fig. 9), the resonance occurs at the lowest field position, 
and the projected distribution of tilts onto the rotation plane is at its broadest, resulting in 
a broader EPR peak. Conservation of the area under the peaks additionally requires that 
the (pE = 0° peak be more intense at its maximum than the <$e = 90° peak. The result is 
that the broad low-field shoulder does not quite get resolved from the main C peak, whereas 
the sharp high-field peak is resolved for several of the resonances — particularly adO and a8. 
In fact, the low-field 4>e = 90° peak contributes significantly to the slow decay of the o;10 
resonance on its low-field side, thus accounting for another previously unexplained feature 
of the Mn 12 -ac single-crystal EPR spectra. We note that recent experiments on deuterated 
samples reveal shoulders on both the low and high field sides of the main resonances.— 

Knowing the relative orientations of the principal axes associated with the E and B\ 
tensors, we can simulate the full hard-plane angle dependence published previously— as 
shown in Figs. 10 to 12. We note that, while the misalignment of the E and B\ tensors has 
previously been suggested from x-ray structure measurements, the EPR studies in ref. l^j ] 
clearly illustrate that the magnetic properties of Mn 12 -ac reflect this misalignment. Fig. 10 
displays the simulated spectra as a function of the field orientation (p relative to one of the 
hard four-fold axes (HBf). These spectra can be compared directly with Fig. 1 of ref. jl^ . 
albeit that the frequencies are slightly different (50 GHz in ref. |l4] and 62 GHz in the present 
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case). The previously published spectra do not include the highest field a 10 resonance due 
to the limitations of the split-pair magnet used in those investigations; one should also be 



careful to note the different labe 
of a8 (highest field peak in ref. 



ing scheme used in this and earlier studies.— The splitting 



14j ) is reproduced very precisely, as are the amplitudes of 



the linewidths (full-width-at-half-maximum — FWHM) variations with angle, as displayed 
in Fig. 11 for peaks a2 to a8. Again, these are the same four linewidths plotted in Fig. 2a 
of ref. Q]. The maximum linewidth or splitting caused by the disorder-induced rhombic 
anisotropy is observed at 4> = ±60° and =f30° relative to the HBf axes (see Fig. 6), and 
the widths are in fair agreement with the £"-value used in the present simulations. The 
orientations of the linewidth maxima, which correspond to the hard and medium two-fold 
axes (HEa and HEb) are, therefore, oriented at ±70° and =f20° relative to the crystal faces. 
The ± and =F symbols are used here to reflect the inversion symmetry of the crystal. We note 
that the linewidth minima are somewhat deeper in the simulation (Fig. 11), as compared 
to experiment.— This is most likely due to other sources of line broadening (D-strain, g- 
strain, etc..) which were not taken into consideration in our model. The main point of this 
simulation is to once again show that the four- fold oscillation in the linewidth/splitting is 
related to the ^-strain, and that the model employed in this study is in complete agreement 
with our earlier investigations. The contour plot in Fig. 12 clearly illustrates the four-fold 
variation in the line positions caused by the intrinsic four-fold B\0\ transverse anisotropy. 
The orientation of one the HBf axes, together with one of the HE axes, is indicated at the 
top of the figure, as is the orientation of the rotation plane corresponding to the data in 
Fig. 5. Overall, Fig. 12 is an excellent reproduction of the contour plot in Fig. 1 of ref. jl^ . 
We remind the reader that the linewidth/shape analysis was employed in ref. [14] to estimate 
the rhombic .E-term associated with the low-symmetry Mni2 variants, and is supported by 
independent magnetization studies.— Contrary to recent assertions,— the EPR linewidths 
are not used in our determination of the easy-axis tilting, which has been inferred entirely 
from the angle-dependent intensities of the EPR spectra. 

Finally, having reproduced the full angle dependence at a single temperature of 15 K, we 
display in Fig. 13 a simulation of the older temperature dependent data shown in Fig. 2. 
We found it necessary to include a small mis-alignment (1°, or 6 = 89°) of the field away 
from the hard plane in order to obtain the appropriate relative intensities of the a and (3 
resonances, i.e. we suspect an approximately 1° mis-alignment in the original experiment 
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(we note that these earlier studies were performed in a cryostat which did not allow for 
in-situ alignment of the sample). Within the hard plane, the simulation assumes that the 
field was applied along the HBf axis (i.e. = 0). Agreement between the simulations and 
the experiment is extremely good. 

V. SUMMARY AND CONCLUSIONS 

Using angle-resolved single-crystal high-frequency EPR measurements, we have clearly 
demonstrated that the molecular easy axes associated with the widely studied Mn 12 -ac 
single-molecule magnet are tilted on a local scale. This tilting provides the much sought- 
after source of the transverse fields necessary to explain the observation of 'odd' resonant 
tunneling steps in hysteresis experiments, i.e. the tilt distribution results in a transverse 
internal field distribution, even when the field is applied precisely parallel to the global 
easy-axis of the crystal. More importantly, on the basis of the angle- dependence of fine 
structures observed in the EPR spectra, we infer that the tilts are constrained along or- 
thogonal planes, i.e. there is a discrete aspect to the tilts, as was previously suspected from 
studies which demonstrated the presence of a locally varying rhombic anisotropy . 14 i 15 This 
behavior is completely reproducible in Mn 12 -ac samples prepared by different methods and 
by different groups . 23 ' 24 Furthermore, our analysis is insensitive to the intricacies associ- 
ated with lineshape analyses, which could be influenced by complex many-body effects, e.g. 
magnetic spin-spin interactions.— Therefore, this study provides compelling support for the 
recent predictions by Cornia et al.,— who have suggested that a discrete form of disorder 
associated with the acetic acid of crystallization provides the dominant source of transverse 
anisotropy responsible for the magnetic quantum tunneling in Mn 12 -ac. Thus, one can rule 
out suggestions that dislocations may be responsible for the symmetry lowering at the local 
scale, since one would not expect the easy axis tilting to be confined to discrete planes under 
such a scenario M 

While our data agree with many qualitative aspects of the Cornia model,— the magni- 
tudes of the easy-axis tilting angles (up to 1.7°), and the rhombic anisotropy associated with 
the tilted molecules (E up to 0.008 cm -1 ), are both significantly larger than the predicted 
values (by a factor of about 5). These results should, therefore, be of significant theoreti- 
cal interest, forming a basis for more refined structure calculations such as those recently 
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published by Park et al.— We tentatively ascribe the observed tilting to the hydrogen- 
bonding between the two acetic acid solvent molecules per Mni2-ac. This conjecture is 
supported by earlier observations that hydrogen bonding must be important to the magne- 
tization dynamics in the Mn 12 cluster because deuteration leads to a significant slowing down 
of the magnetization reversal dynamics as measured by high-frequency ac susceptibility.— 
We note also that all of the unusual solvent-induced anomalies reported in this paper for 
Mni2-ac are absent in the EPR spectra for a related high-symmetry (S4) Mni2 complex 
[Mn 12 Oi 2 (0 2 CCH 2 Br) 16 (H 2 0) 4 4CH 2 Cl 2 ] which possesses a full compliment of four CH 2 C1 2 
solvent molecules per Mni 2 , and none of the low symmetry hydrogen bonding environments 
present in Cornia's model. 26 
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Figure captions 

FIG. 1. (a) Energy level diagram for Mrii2-ac with the field applied in the hard plane 
{9 = 90°, 4> = 45°); the crystal-field parameters used for the simulation are given in the main 
text. The labels on the left hand side of the figure correspond to the spin projection along z 
(low- field m z basis), while the labels on the right correspond to the spin projection along the 
applied field direction (high-field m x basis). The vertical blue sticks illustrate the origin of 
the a EPR transitions, while the filled red circles represent the high-frequency ,5-transitions. 
(b) A comparison between measured and calculated a and /3 EPR transition frequencies, as 
a function of the applied magnetic field strength. The open blue squares (a) and red circles 
(/3) represent experimental data obtained at 8 K. The thick curves illustrate the expected 
behavior of the EPR transition frequencies for precise field alignment within the hard plane 
(9 = 90°), i.e. no tilting of the molecules. Note that, for this case, the a frequencies decrease 
smoothly to zero as the field tends to zero, whereas the (3 frequencies go through a minimum 
(~ 90 — 95 Ghz for the highest field transitions) and then tend to a finite zero-field offset. 
The thin curves illustrate the dependence of the a and f3 transition frequencies upon tilting 
the field away from the hard plane; each thin curve represents an additional 0.5° of tilt away 
from the hard plane (9 = 90°), and the curves have been truncated at low fields in order to 
simplify the figure. The red arrow illustrates the depression of the minimum f3 frequency, 
while the blue arrow illustrates the corresponding increase in the minimum a frequency. 

FIG. 2. Temperature and frequency dependence of the EPR spectra obtained with the 
magnetic field applied approximately parallel to the hard four-fold axis of the crystal ; 30 ! 31 the 
frequencies and temperatures are indicated in the figure. These experiments were performed 
using instrumentation that did not allow for in-situ alignment of the sample. Based on 
subsequent analysis, we estimate a slight (~ 1°) mis-alignment of the field away from the 
hard plane (see also Fig. 13), i.e. 9 = 89° and <fi ~ 0°. The data clearly demonstrate that 
both the a and /3 resonances persist to the lowest frequency (44.3 GHz), in apparent conflict 
with the predictions of Fig. 1, thus suggesting that a significant fraction of the molecules 
may be tilted with respect to the global easy axis of the crystal (see text). The 7 resonance 
has been previously attributed to a small fraction (< 5%) of Mni2 molecules having lower 
symmetry (i.e. different crystal field parameters) than the acetate.^ - 
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FIG. 3. Simulations of 15 K, 44 GHz, EPR spectra as a function of the field orientation 9 
(in 0.2° increments) relative to the easy axis of a single molecule ((f) = 0°). In (a) the vertical 
scale represents absorption (arb. units) while, in (b), absorption is indicated by the darker 
shaded regions. The parameters used for these simulations are given in the main text. In 
(b), the dashed curves represent constant values of the longitudinal field Bn. The key point 
to note is the alternate appearance and disappearance of even a peaks and odd (5 peaks. 

FIG. 4. Simulations of 15 K, 62 GHz, EPR spectra as a function of the field orientation 9 
(in 0.2° increments) relative to the easy axis of a single molecule ((f) = 0°). In (a) the vertical 
scale represents absorption (arb. units) while, in (b), absorption is indicated by the darker 
shaded regions. The parameters used for these simulations are given in the main text. In 
(b), the dashed curves represent constant values of the longitudinal field B\\. The key point 
to note is the alternate appearance and disappearance of even a peaks and odd (3 peaks. 

FIG. 5. a) Experimental EPR spectra obtained at 15 K and 62 GHz, as a function of 
the field orientation 9 (in 0.18° increments) relative to the global easy (z-) axis of a single 
crystal sample; the red data represent 9 = 90°. These experiments were performed using 
instrumentation that did allow for in-situ rotation of the sample about 1 axis. 39 We estimate 
(see main text) that the rotation plane was 10° away from one of the hard B\ axes within 
the hard (xy-) plane. To aid comparisons, the data in b) are plotted in the same fashion as 
Fig. 4a, while the data in c) are plotted in the same fashion as Fig. 4b. In (b) the vertical 
scale represents absorption (arb. units). In (c) absorption is indicated by the darker shaded 
regions; 15% and 50% (of maximum absorption) contour lines have been included to aid the 
discussion in the main text. No noise reduction/filtering has been applied to the data. The 
main points to note are the significant overlap of intensity due to the a and (3 resonances, 
and the distinct shoulders on the high field sides of the peaks which have been emphasized 
with the 50% contour lines. 

FIG. 6. Schematic illustrating the discrete nature of the proposed tilting of various Mn 12 
species. Upper panel: illustration of the tilt distribution for the A molecules, wherein 
the tilts are confined to the plane containing the crystallographic easy (z-) axis and the 
hard two-fold axis of the B molecules {HEb = medium two-fold axis of the A molecules); 
meanwhile, the hard two-fold axis of the A molecules (HE a) is perpendicular to the tilt 
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plane; an orthogonal situation exists for the B molecules. Lower panel: a polar plot (9 is 
radial and (f> is azimuthal) illustrating the orientations of the two tilt planes (shaded areas), 
and the principal axes associated with the intrinsic four-fold anisotropy (Of) and disorder 
induced two-fold anisotropies (Of) associated with two of the Mn 12 variants (A and B); 
H Ea and HEb refer respectively to the hard two-fold axes for the A and B molecules; HBf 
and MBf refer respectively to the global hard and medium four-fold axes; the purple arrow 
indicates the plane of rotation for the experiments shown in Fig. 5, while the red and blue 
dashed lines respectively indicate the extent of the projections of the tilts of the A and B 
molecules onto the rotation plane. Note the = 30° offset between the Of and 0\ tensors. 

FIG. 7. Simulated contributions to the angle dependent EPR spectrum for the A, B 
and C molecules, for rotations in a plane which is inclined A0 = 20° away from the tilt 
plane associated with the B molecules (70° away from the A tilt plane, see also Fig. 6). The 
temperature is 15 K and the frequency is 62 GHz; see main text for a detailed explanation 
of these simulations. The E'-terni associated with the A and B molecules splits their contri- 
bution to the spectrum. The near coincidence of the rotation plane and the hard axis of the 
A molecules shifts the contribution to the intensity to higher (lower) fields for the A (B) 
molecules. Meanwhile, the near coincidence of the rotation plane and the B tilt plane (see 
Fig. 6) projects almost the full tilt distribution for the B molecules, i.e. the absorption due 
to the B molecules extends over a wider angle range either side of the hard plane (9 = 90°) 
compared to either A or C, thus accounting for the overlapping a and (3 resonances in 
Fig. 5. In contrast, the angle dependence of the EPR spectra associated with the A and C 
molecules does not differ greatly from the simulation in Fig. 4 (with no tilting). 

FIG. 8. a) Angle dependent 25%:25%:50% summations of the normalized spectra for the 
A, B and C molecules from Fig. 7; each successive trace corresponds to a 0.2° increment in 
9, and the red trace represents 9 = 90°. In b), the same data are displayed in a 3D gray- 
scale for direct comparisons with earlier figures; the 10% and 50% (of maximum absorption) 
contours are intended as guides to the eye. These simulations are to be compared with the 
experimental data in Fig. 5. In particular, the high-field shoulder is reproduced (see also 
Fig. 10), as well as the considerable angle overlap of the a and (3 resonances. 

FIG. 9. Simulated evolution of absorption due to either of the tilted species (A or B) as 
a function of the field orientation (cf) E , within the hard plane) relative to the hard two-fold 
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axis for that species (note - this is not the same as the angle <p); the temperature is 15 K 
and the frequency is 62 GHz. With the field along the hard two-fold axis HE = 0°), the 
resonance occurs at the highest field position, while the projected distribution of tilts onto 
the rotation plane is at its narrowest. Consequently, the EPR peak is narrow, forming a sharp 
shoulder on the high field side of the summed spectra (see Figs. 8 and 10). Conversely, with 
the field along the medium two-fold axis {<pE = 90°), the resonance occurs at the lowest field 
position, and the projected distribution of tilts onto the rotation plane is at its broadest, 
resulting in a broader EPR peak, which contributes to the low-field tail of the summed 
spectra. 

FIG. 10. Simulated ^dependence of the 62 GHz, 15 K, EPR spectra incorporating 
the A(p = 30° mis-alignment between the intrinsic four-fold 0\ tensor and the disorder- 
induced two-fold 0\ tensors for the A and B molecules. These data compare favorably with 
previously published data (see Fig. 1 of ref. [l^j]). 

FIG. 11. Azimuthal angle ((ft-) dependence of the full-widths-at-half-maximum (FWHM) 
for several of the a peaks in Fig. 10 (ct2 to ct8). This figure is to be compared directly with 



Fig. 2a of ref. 



FIG. 12. Grayscale contour plot corresponding to the data in Fig. 11. The black dashed 
line indicates the orientation corresponding to the data in Fig. 5. The orientations of one of 
the HE axes and one of the HB\ axes are indicated at the top of the figure. This figure is 
to be compared directly with the contour plot Fig. 1 of ref. 14] (note the slightly different 
labeling used in these two figures^) ; the white dashed lines in the present figure correspond 
to the orientations of the two experimental curves in Fig. 1 of ref. jl^j ]. 

FIG. 13. Full simulations of the temperature and frequency dependent experimental data 
in Fig. 2 with the magnetic field applied approximately parallel to the hard four-fold axis of 
the crystal; the frequencies and temperatures are indicated in the figure. In order to obtain 
the best agreement, a slight 1° mis-alignment away from the hard plane was taken into 
account, i.e. 9 = 89° and (p ~ 0°. Indeed, the relative intensities of the a and f3 resonances 
are extremely sensitive to the angle 6. This is easily understood upon inspection of Fig. la, 
where it is seen that the relative energy separations of the lowest lying states are extremely 
sensitive to the field orientation. Therefore, the corresponding transition probabilities, which 
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depend on the differences between Boltzmann factors, will also depend very sensitively on 
the angle 9. 
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FIG. 1: S. Hill et al. 
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FIG. 2: S. Hill et al. 



28 




FIG. 3: S. Hill et al. 
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FIG. 4: S. Hill et al. 
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FIG. 5: S. Hill et al. 
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FIG. 6: S. Hill et al. 
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FIG. 7: S. Hill et al. 
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FIG. 8: S. Hill et al. 
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FIG. 9: S. Hill et al. 
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FIG. 10: S. Hill et al. 
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FIG. 11: S. Hill et al. 
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FIG. 12: S. Hill et al. 
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FIG. 13: S. Hill et al. 
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